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Abstract 

We formulate the quark sea contributions to the energy and various physical 
observables in terms of Euclidean Green functions and their K-spin partial 
wave reduction. In this framework it is not necessary to discretize the con- 
tinuous spectrum by introducing a finite boundary. Using this formulation we 
perform a new self-consistent computation of the nucleon state in the Nambu- 
Jona-Lasinio model. Besides this technical advantage such an alternative com- 
putation scheme makes it possible to obtain the numerical predictions of the 
model in an entirely independent way. We use the Pauli-Villars cutoff to de- 
fine the model. Our results for the nucleon energy, the mesonic profile, and 
various observables essentially confirm results obtained by other groups using 
this regularization. The results for g\, obtained on the basis of quark currents, 
are close to the experimental value. 
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1 Introduction 



The Nambu-Jona-Lasinio (NJL) model |J has received a wide attention dur- 
ing the last decade [f| , following previous attempts || |J to describe the nu- 
cleon in the large- N c limit on the basis of mesonic degrees of freedom. In the 
NJL model, the mesons appear as effective degrees of freedom, parametrizing 
condensates of the basic fermion fields. The basic model is a quark model 
with a four-fermion interaction, and therefore nonrenormalizable. The ultra- 
violet divergences are handled by introducing a cutoff which stays finite. Its 
functional form and numerical value, therefore, are relevant for the predic- 
tions of the model. Various cutoffs have been introduced, a good review of 
the different possibilities and their numerical impacts is given in ||. Previ- 
ously, the Schwinger proper-time cutoff was used almost exclusively until it 
was found that this regularization violates the momentum sum rule for 
the parton distributions, and that Pauli-Villars regularization is favored in 
this respect. 

The use of the Pauli-Villars regularization opens the possibility to in- 
troduce a numerical technique for computing the effective action which has 
been developed previously [0, § and has been applied to a various physi- 
cal problems involving fluctuation determinants [[], [ITJ. It is based on using 
the Euclidean Green functions instead of summation of levels. Computing 
functional determinants and other expectation values involving the quark 
continuum by summing over levels and using Minkowski space wave func- 
tions requires the introduction of space boundaries in order to discretize the 
spectrum. The associated space cutoff has to be removed, introducing a 
numerical limiting procedure. This limiting procedure seems to be techni- 
cally well under control, as can be seen, e.g., by comparing computations of 
the sphaleron determinant using level summation ]TT[ and using Euclidean 
Green functions M. Nevertheless, level summation is certainly not a very 
economical technique. At the same time an alternative computational ap- 
proach presents the possibility to obtain the predictions of the model in an 
independent way. 

The idea of replacing level summation by integrals over Euclidean Green 
functions is actually a rather old one. It seems to go back to Wichmann 
and Kroll [12[]. It has been used extensively in calculations of the vacuum 
polarization in strong fields [|1^, [13j and of the Casimir effect [14]] . A related 
way of encompassing the discretization of the spectrum has been proposed 
by Moussallam [JTH], who uses Minkowski space phase shifts. 
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There are various methods to perform self-consistent computations |TB 



As the equation of motion for the meson profiles requires to solve the equa- 
tion of motion for the meson field, or, equivalently, the chiral angle, it is 
advantageous to be able to compute the functional derivative of the effec- 
tive action with respect to the meson field. A technique for computing such 
derivatives using Euclidean Green functions has been set up recently [17] and 
used for a self-consistent computation of the bubble nucleation rate in the 
electroweak theory ]18[ . It is the purpose of this work to transfer these tech- 
niques, mutatis mutandis, to the NJL model. Apart from the self-consistent 
computation of the mesonic profiles of the nucleon we also formulate the 
sea quark contributions to various other observables, such as the moment 
of inertia, in terms of Euclidean Green functions. This latter aspect of our 
work is of course important; without it, one would have to go back to level 
summation in computing these observables and our technique would loose 
much of its attractiveness. 

The NJL model, as considered here, is not a unique theory; there are 
many versions of it, the most elaborate ones |19[ include the p, a\ and u 



vector mesons as well, as it was done before in the Skyrme model [^0| |21 
However, even within the restricted class of models in which only the n — a 
fields are taken into account, there is a wide variety. These models differ by 
the kind of regularization, as mentioned above, but also by applying it to 
various parts of the spectrum. In a renormalizable theory one would apply it 
only to the lowest perturbative contributions; in the NJL model it is applied 
to the finite parts as well. Furthermore, it was usually only applied to the 
quark sea , and not to the valence contribution. Some models also differ by 
the way in which the meson degrees of freedom are varied: the variation can 
be extended to both a and tt fields, or be restricted to the "chiral circle", 
the 4-sphere being defined by a 2 + 7r 2 = f 2 . Unfortunately, the existence 
of nucleon solutions is not a robust property of the theory, such solutions 
are found only within a small subset of the different versions, a situation 
that must be considered as unsatisfactory. We will not add to the ongoing 
discussion (see e.g. P^|), we will use that version of the model recently used 
by Pobylitsa et al. [23| for computing the parton distributions. It only takes 
into account the 7r and a fields, which are varied on the chiral circle only, 
and Paul-Villars regularization is applied to the quark sea, not to the filled 
bound state. 
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2 The model 



Starting with the NJL-Lagrangian [ITJ 



C 



G 



one obtains after the standard bosonization procedure 



Si 



NJL 



i'y fl d^-g(cr + in-T>y 5 



2\ . 



(2.1) 



(2.2) 



The parameters of the model are the fermion self-coupling G, the quark mass 
m, and the cutoff scale A. In the bosonized version these parameters appear 
as the quark-meson coupling g and the symmetry breaking mass parameter 
ji. They are related to the basic parameters as g = fiyG and mfi 2 = gf^m 2 . 
The latter equation expresses \i in terms of physical constants and of the 
coupling g which remains a free parameter. A further relation is obtained 
from the gradient expansion of the effective action. The resulting kinetic 
term of the pion field is normalized correctly if the Pauli-Villars cutoff is 
fixed as 



A = M 



\ 



exp 



471-2 



(2.3) 



If the a — 7r field is varied only on the chiral circle, the second term in Eq. 
( |2.2j ) is absent. For static pion fields the action is proportional to the time r; 
the two remaining parts of the effective action then contribute to the energy 
as 



-Efer 
El 



1, 



Tr log -i"/^ + M(x) 



br = -mlf n J d 3 x (a - f n ) . 



(2.4) 
(2.5) 



Here the trace is taken over the quark sea and - in the case of baryons - over 
the filled bound states. M(x) is given - on the chiral circle - by 



M(x) = g a + rysT ■ 7r 

= M exp {«75T • <p(x)} 
= M cos(\4>(x)\) + ry 5 r ■ <f>(x) sin(|0(a;)|) 



(2.6) 
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where we have introduced the "dynamical quark mass" M = gf n . With the 
hedgehog ansatz 4>(x) = xd{r) the mass becomes 



Mix) 



M 



cos(i9(r)) + ij$T ■ x sin(i?(r)) 



(2.7) 



3 Basic relations 

Given the profile i?(r), the energy of the corresponding nucleon state con- 
sists of the symmetry breaking part that can be evaluated trivially, and the 
contributions of valence and sea quarks. In order to evaluate the valence 
quark contribution we have to find the bound state energy by solving, with 
appropriate boundary conditions, the Dirac equation 

(iu - H)^ (x) = . (3.1) 

Here we have introduced the Dirac Hamiltonian 

H = -ia ■ V + -foM(x) . (3.2) 

The computation of the quark sea contribution is more involved. We will 
recall here a method introduced previously [[jj in which the computation of 
the zero point energy is related to the Euclidean Green function 

Se(x ,* >) = £«?MM (3.3) 
V -w + E a 

The subscript a is a formal notation for the discrete and continuum eigen- 
states of the Dirac Hamiltonian H; we also indicate the positive energy eigen- 
states by a > 0, the negative ones by a < 0, and the valence eigenstate with 
a = 0. Se satisfies the equation 

(iu - H)Se(x, x r , v) = -5 3 (x - x 1 ) . (3.4) 

The zero point energy 

E sea , = ^2 (3.5) 

can be computed as a contour integral around the positive imaginary axis in 
the complex //-plane, see Fig. 1, as 

£sea= / ^UTI [ d 3 xS E (x,X,ls) . (3.6) 

JC- 2m J 
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Deforming the contour to run along the real v axis, and subtracting the 
zero point energy of the free Dirac operator H = — ia- V + 70M, the integral 
takes the form 

/OO Ajj r 
^— z/Tr j d 3 x[S E (x,x,u) - S Efl (x,x,u)] . (3.7) 

The original contour, and therefore also the deformed one, takes into account 
all negative energy states, whether they are continuum or bound states. The 
only state that has to be considered separately is the positive energy bound 
state that is occupied by the valence quarks and, therefore, should be included 
as well. This can be done by deforming the contour to the one presented 
in Fig. 1 as a dashed line. It is, however, more convenient to write the 
contibutions of this state separately, as we will do here. 

It is convenient to introduce the bosonic Green function G E via 

S E = {%v + H)G E (3.8) 

which satisfies 

[v 2 + H 2 ) G E (x, x', u)=[u 2 -A + M 2 + V(x)] G E (x, x', v) = 5 3 (x - x') 

(3.9) 

with the potential or vertex operator 

V(aj) = iy ■ VM(x) . (3.10) 
In terms of G E the energy can be written as 

E = / — v 2 \ ciVTr [G E (x, x, v) - G Efi (x, x, u)\ . (3.11) 

JO IT J 

The expressions for the zero point energy and the subsequent manipulations 
are formal. Even after subtracting the free zero point energy, they only make 
sense if properly regularized. In order to understand the divergences of the 
subtracted zero point energy we use the resolvent expansion of the Green 
function with respect to the potential V; if it is inserted into the expression 
for the zero point energy we obtain 

E o = E hr U d 3 ^G E ,o(x t - u)V(xi) (3.12) 

n=\ 2n i=l 
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with x n = xq. We have subtracted the free zero point energy by omitting the 
zeroth order in the sum over n. The first order term vanishes after taking 
the trace, and the only divergent term is the second order one. Explicitly, it 
takes the form 

4 J (2it) 6 Jo J (27r) 4 [p 2 + M 2 + q 2 x(l - x)\ 

(3.13) 

where V(q) is the Fourier transform of the potential. The logarithmically 
divergent integral can be defined by using, below the integrand, the Pauli- 
Villars subtraction 

^ = I / JffiS^M-ti Jo dx Tib i lnA2 ~ ln K + q2%{l - X) l } • 

(3.14) 

Unlike in the case of renormalized perturbation theory here the regularization 
will be extended over the finite contributions as well, so that the regularized 
zero point energy reads 

£ ,rcg = E {M) - ^E (A) , (3.15) 

where the subtraction is to be understood to be done below the x and v 
integrals, respectively (see below), before the partial wave summation and v 
integration. The factor M 2 /A 2 takes into account that the potential V(g) 
contains a factor M, and a factor A in the subtracted part. These prefactors 
have to be compensated in order to ensure the cancellation of the divergent 
integrals. 

The perturbative expansion, Eq. ( p.!2j ), can be used JL7| to obtain an 
expression for the derivative of the zero point energy 



= y(~l) n fdu r 3 5V{x x ) 

Ha{z) tl 2 J 27T J Xl 50 a (z) 

/ (il d 3 XiGE,o(xi - Xi^i, v)V{xi)^ G E ,o(^i - x n , v) . 

(3.16) 



\i=2 



The perturbative sum can be recollected into the full nonperturbative Green 
function, so that 



- tr W d 3 x^- f°° -G~\x x u) (3 17) 
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Here and in the following we use superscripts (j) to indicate that the ex- 
pression is of order j in the potential V, or its derivative. The symbol (j) 
indicates that the expression is evaluated to all orders in the potential start- 
ing with order j. So the exact zero point energy, after subtracting the zeroth 
order, and taking account of the vanishing of the first order, is of order (2), 

as indicated on the l.h.s. of the last equation. on the r.h.s. includes all 
orders but the zeroth one. It remains to evaluate the derivative of V(x) with 
respect to <p a . One obtains 



6V(x) 



iM~f ■ VJ 3 (x - 
+il5 (r a - (j) a (x) t ■ 4>{x 



iy 5 r a exp(iy 5 T • <j){x)) (3.18) 
sin \ 4>(x) 



Inserting the hedgehog ansatz, <p(x) = ah?(r) and using the fact that within 
the trace the expectation value of r must be parallel to 0, i.e., r a — > x a T ■ x, 
the derivative of the zero point energy takes the form 



6E, 



(2) 



-Mtr7 5 r a cos(i?(r)) — iy$T ■ zsin({}(r)) 



4>(X)=Zd(r) 



dv 



Z7T 



(3.19) 



The gradient acts on the Green function at equal arguments. It is taken after 
carrying out that limit. The Euclidean Green function can be expanded f24 
with respect to K-spin harmonics ^n' Kz (f° r details see Appendix A) as 



G E (z,z' v) 



E 9- 

K,K Z ,P 



^ p (r,r>,v)Z«> K >(z)®Z^(z') 



(3.20) 



The radial Green functions form 4x4 matrices; they can be written in terms 
of mode functions^] f^ + (u,r) and f%~(v,r) which are solutions regular at 
r = and as r — > oo, respectively, of a system of radial differential equations 
given in Appendix A. Explicitly, they are given by 

g mn (r,r', v) = « [0(r - r>)fZ + (is,r)f:-(v,r>) + 6{r' - r)f^(u } r)f: + (u,r') . 

(3.21) 



In the following we omit the K-spin and parity superscripts. 
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The superscript a labels 4 linearly independent solutions. 

In this basis, and using the reduced Green functions, the zero point energy 
takes the form H 



E, 



(2) 



N c . 



dv 



7T 



— v 2 \ drr 2 J2(2K + 1) 



K,P 



(2) , (2) , (2) . (2) 
011 + 922 + 933 + 9il 



(3.22) 



where the Green functions are taken at r = r' . Analogously, the functional 
derivative of the energy is obtained as 



SEP 



(p(Z)=Z$(r) K,P 



x / c?v< sia(i?(r)) 



(2K + 1; 



— cos(i?(r)) 



1 ( (i)' , (i)' (i)' , (i)' 

2\~9U + 922 - 933 + 044 



+^[(K- l)g$ + (K + l)g% + Kg&> + (K + 2)g M 



(i) 



+2J^+1)(-^ - 



3^ 



(3.23) 



Both expressions have to be regulated as implied by Eq. ( |3.15 ). 

The NJL-soliton is a system with baryon number equal to one. Therefore, 
one has to add the bound state part of the fermionic energy 



^comp = NcE bou 



The eigenvalue equation for the bound state reads 



-A + V 0+ (a;) 



(3.24) 



(3.25) 



For K p = 0+ the spinor ipQ is determined by two radial wave functions ho(r) 
and j'o(t), corresponding to the components M3 and M4 for K p = + and 
for the bound state energy Eq; the potential V 0+ is a 2 x 2 matrix given 
in Appendix A. The eigenfunctions are normalized as / \ipQ(x)\ 2 d z x = 1. 
Differentiating the bound state equation with respect to 4> a (z) and projecting 
with ipQ one finds 
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-LOq 



2uq 



X 



(3.26) 



Inserting the explicit expressions for the potential and the eigenfunctions, 
the derivative of the bound state energy takes the form 



5E ho 



M- 



4>(Z)=Z'&(r) 
N, 



z a < sin($(r)) 



K( r )Jo(r) +j'o(r)ho(r) + -h (r)j (r) 



cos(i?(r)) 



-h' (r)h (r) +f (r)j (r) + -f (r) 



. (3.27) 



The bound state energy is convergent. Thus it does not need to be regular- 
ized. With a finite regulator, however, this is subject to some arbitrariness. 
The same argument would hold for any finite subset of sea quark states, or 
of all the finite parts higher order terms of the perturbative expansion. 

Finally, the mesonic part of the energy and its derivative have to be 
evaluated. This is straightforward, as they are given by simple analytical 
expressions. On the chiral circle one finds 



E 



br 



—Anm^f^ drr cos(i?(r)) 



and 



5E 



br 



(f)(Z)=Z'd(r) 



m 7 ,f 7r z a sia('&(r)) 



(3.28) 



(3.29) 



4 Perturbative expansion of the Green func- 
tion 



After reduction of Eq. ( |3.9| ) to K-spin partial waves (see Appendix A) the dif- 
ferential equation for the partial wave Green functions g m n{ r , r ', v) becomes 



nk \ dr 2 ^ r dr 
~- -8nm<5(r - r')/r 2 



2 d K n {K n + 1) 



k 2 - V nk {r) 



gkm(r, r', v) 



(4.1) 



where k = \/v 2 + M 2 and where, again, we suppress the partial wave indices 
K and P. The potential V mn depends on the K-spin, its explicit form is 
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given in Appendix A. As already mentioned, we use for the numerical com- 
putation of the Green functions their standard expression (|3.21| ) in terms of 
mode functions f 



The functions form 4x2 linearly independent 
systems (index a±) of 4-component solutions (subscript n). A form indepen- 
dent of the choice of basis is given in f|, H3 • Here we use a special convenient 
basis; it is defined by splitting off the free solutions, i.e. the modified Bessel 
functions b~^ n (nr) = /cx n (^ r ) and b^ n (nr) = ix n {^r) via [] 



and by imposing the boundary condition 

lim hZ^iv, r) 







(4.2) 
(4.3) 



The truncated mode-functions are obtained by solving the equations 
- (I b r± 

dr 2 



KT 



— + 2 ( - + ^ttt— 

bis (nr 



d_ 

dr 



V K (r) 



a±. 



S a + h 



z/,r) 



b% 



b Kn M 



(4.4) 

or, in a short form, 

Dh = V(l + h) . (4.5) 

This equation can be used for a perturbative expansion. Obviously, the 
functions h^{v, r) vanish to zeroth order in V, so, in the notation introduced 
in the previous section, they are of order (1). Once these solutions are known, 
the differential equation may be iterated to obtain the contribution of order 
(2) via 

DhP = V (l + hP) , (4.6) 
Dh& = Vh®. (4.7) 
In terms of these functions the expression (J3.21) becomes, for r > r', 



K 



h£ )m -(v,r') + h$ n+ (v,r) + h^ a ~(u,r')h^ a+ (u,r) 



t( 2 ) 



xk 

K 



K m (Kr i Kn (Kr 



(4.8) 



xk 



K m (Kr i Kn (Kr 



(4.9) 



4 As for the mode functions /" we omit the K - spin and parity assignment for the 
truncated mode functions h"^. 
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These expressions are ready for being inserted into Eqs. ( p.22|) and ( |3.23| ). In 
renormalized perturbation theory one would, using further iterations, reduce 
these expressions to order (3) and evaluate the second order analytically as 
in Eq. ( gig) . 



5 Quantization of collective coordinates 

The hedgehog system is invariant under K-spin, i.e., under combined space 
and isospin rotations. For a rotating hedgehog state the Hamiltonian is 



modified as 25 



tf n (0) = H(x) - -n • r = -iot ■ V + -f M(x) - hi ■ r , 



(5.1) 



where ft is the angular velocity. Assuming O to be small, the problem can 
be treated perturbatively. The first order in O vanishes, in second order one 
obtains 

E (n) + ln a 9 ab n b ] , (5.2) 



S e f[ = —r 
where Q n b is the moment of inertia 



e, 



ab 



(5.3) 

= 95 a b- The collective coordinate 



9 a b is proportional to the unit matrix, 9 a b 
can be quantized in the usual way, leading to an extra term J(J + l)/29 in 
the energy. 

Taking the second derivative with respect to Q of the fermion Green 
function 



SS(x,xf,v) = -lx\ 



iv-HnW, 



\x') 



(5.4) 



one obtains 

6 2 S%{x,x',v) 



—\x\ 



1 



1 



\X\ 



1 iv - H n (fl) 2 iv - iZ„(n) 2iv- H n (Q] 
iv - #n(n) 2 iv - fln(O) 2iu- H n (fl) 

x'). 



\x') 
\x') 



id 1 1 

4d^^ Ta tv - H n (fl) Tb tv - H n (Q) 



(5.5) 
(5.6) 
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Inserting this equation in Eq. (3.11) the tensor can be calculated via 



17=0 



dv 



57T 



d 3 x / d 3 x'tr r a S E (x, x', v)r b S E (x', x, v) 



(5.7) 

where S E (x,x', v) is defined by Eq. ( |3.4| ). This expression can be rewritten 
with the "bosonic" Green function (I3 - 



e. 



ab 



dv 



S7T 



xtr 



d 3 x J d 3 x' (5.8) 
r a {iv + H(x)} G E {x, x', v)r b {iv + H(x')} G E {x' , x, v) 



The bound state is occupied and, therefore, included into the negative con- 
tinuum by choosing the dashed integration contour displayed in Fig. 1. If 
the Green function is expanded into eigenfunctions of the Hamiltonian as 



in Eq. (IP), it can be readily verified that this way one obtains transitions 
between the positive and negative continuum states as well as transitions be- 



tween the bound state and the positive continuum as discussed in |26| . Here 
this detailed structure is not explicit. 

As for the energy, we will consider the bound state contribution sepa- 
rately. We decompose the dashed contour into a contour running along the 
real v axis and a small circle around the bound state pole. We will denote the 
former contribution which describes the continuum-continuum-transitions by 
a superscript c — c, and the bound state contributions by the superscript b — c 
as it involves matrix elements between the bound state and the continuum. 
We begin with considering the continuum contributions. 

The partial wave reduction of the Green function Ge{x,x',v) has been 
introduced above and is discussed in Appendix A. Here we need the partial 
wave reduction for the fermion propagator = {iv + H)G% we denote 
the associated matrix elements by s^ n {r,r' ,v). This Green function can 



again be decomposed into mode functions as in Eq. |3.21| , one just has to 
replace the functions f n {u, r) by the fermionic mode functions u n {v, r). These 
are given explicitly in Appendix A. Furthermore, we have to take the trace 
with isospin matrices r a . As 9 ab oc 5 a b it is sufficient to compute 9 33 . The 
action of r 3 on the K-spin harmonics is given in Appendix B. Using these 
expressions, the angular integration over dQ r and dQ r r and the summation 
over the third component of K and K' can be performed. K' is fixed by the 
angular integration to the values K, K — 1 or K + 1. Finally, one obtains 



dVt r / dClr'tr 



r 3 {iv + H{x)} G E {x, x', v)t 3 {iv + H{x')} G E {x', x, v) 
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E 

P,K=1 



{K + 1){2K + 1) 
3K 



s n (r, r',u)s u (r',r, v) + s 22 (r, r' ,v)s 22 (r' ,r, v) 



+s 12 (r, r', u)s 21 (r', r, v) + s 21 (r, r', u)s 12 (r', r, v) 



(2K+l)K( Ki 



V- 11 I K / I \ Ki I \ i Ki I \ Ki I \ 

3{K +T)~ r 33 ^' r ' ^ S33 ^ r ' r ' ^ + s 44l^ r , fJ^O , r, v) 

+4( r > r', i/)s^(r', r, z/) + s^(r, r', u)s^(r', r, z/) 

(2 ^ 3 +1) {4(^ r', u)s^(r', r, v) + 4(r, r', i/)s£(r', r, v) 

+sf 3 (r } r',u)sf 1 (r f ,r, v) + 4( r > r ' , u)s^(r' ,r, u) 
+s*(r, r', z/)4(V, r, i/) + s^(r, r', ^)4( r ' ; r, i/) 



4K 2 - 1 



'32V 5 ' 5 "J°23^ 5 ' 5 

+4( r ? r', ^)si(r', r, u) + 4( r > r', i/)s^(r', r, i/) 

4( r > r', z^)4f r, i/) + sg(r, r', v)s« -1 (r', r, z/) 
+4(r, r', z/)4- 1 (r', r, i/) + s£(r, r', z/)4~ 1 (r', r, i/) 

+ s 33 _1 ( r 5 r 'i y ) S ll ( r ',r, U) + 4~ 1 ( r ) r ' ) I/ ) s 22( r ' ; r ) V) 



+4" 1 (r, r', r, v) + 4" 1 (r, r', z/)4(r', r, i/) 



(5.9) 



As we need this expression in order (2), the modified Green functions have 
to be inserted in order (1). One can take advantage of the factorization of 
the Green functions into mode functions, see |A.17j , to rewrite this expression 
in the form 



C 



S 2 E 



°° dv ~ 
-oo 47r pK=l 



drr / dr r 



I '2 



H 



IK ^ 



za r 



2K 



(if + l)(2if + 1) 
3K 

j K(2K + l) Hrf3a 



ik \ 



2K + 1 



2K + 1 



2K \"i 



sk y^i 



UK 2 -I 



za r — 



-/3a, 



z/, r 
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v, r 
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(z/,r) 



AK 2 - 1 



3K 



The functions H^ 1 (i/,r) are defined as 



tt±<xP, 

I7±<xP, 
U 2K 1 



2/,r) 
z/,r) 



u,r) 



n AK 



r)/££(i/, r) + «§*(„, r)/££(i/, r) 



f/3± 



■it 



,a± 



(5.10) 



(5.11) 
(5.12) 



"S^.O/S-x^O +«S(«/,r)/g_ 1 ( I /,r) ,(5.13) 



«± 



i(^r)/S(^r) + <±_ 1 (u,r)fi%(u,r) ,(5.14) 



0±, 



in terms of the mode functions /" ± and w~ defined in Appendix A. One has 
to combine the orders of H^ af3 in such a way that the result is of total order 

(2): 



flC-c(l) 



H+UH-W + ff+Mtf-d) + H HD H -(o) . 



(5.15) 



In fact the first order part which has been included on the right hand side 
for practical convenience vanishes. The functions if are obtained as 



jy±(l) _ J±(l) J±(l) + yr±(0) J±(l) + J±(l) y 



±(0) . 



(5.16) 



the functions H ±( - ^ are composed of free Bessel functions. 

Since the imaginary part of the integral of Eq. ( |5.10| ) is antisymmetric, the 
result equals twice the real part, integrated from v = up to oo. It is implied 
that the expressions have to be regularized by Pauli-Villars subtractions. 

Having presented the continuum-continuum contributions to the moment 
of inertia we now turn to the bound state contribution. This contribution is 
given, in terms of eigenf unctions of the Dirac operator by [^6|, |27[ 



d 2 u 



dtt a d{l a 



Nc 
2 



m^bou 



JTJrn _ JTJbou 



(5.17) 



Using the (|3.3|) for the Green function we find that this expression for the 
moment of inertia is identical to 



nb-c 
U ab 



8 2 UJ 



dVl a d{l 



2 



tr J d 3 xd 3 x'M x )^o( x 'KS(x',x,-tE hou )T b . 

(5.18) 
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The Euclidean Green function at any imaginary argument can again be re- 
lated to the bosonic Green function via 



S E (x, x', -iE hou ) = {H + E bou )G E (x, x', iE 



bou\ 



->bou\ 



(5.19) 



Since the energy of the bound state is smaller than the mass M, the calcu- 
lation of the Green function is analogous to the one for the continuum part, 
Eq. (ETT0D , with 

n 2 = M 2 -E 2 hou . (5.20) 

Furthermore, the valence state + only couples to K p = 1 + continuum states. 
Therefore, the expression for this contribution reduces to 



e. 



b-c 
ab 



N r 



poo 




\ drr 2 


nt a M f 


Jo 


Jo 



[ dr'r ,2 m a (p, r 1 ) 
Jo 



+ni a [ 



with 



Hf a {v,r) 



K 



dr'r' 2 Hz a {v, r' 



ufflv, r)h {r) + uffiv, r)j {r) 



H± a {v,r) 



2-E'bou 



K 



m(v,r)h (r) + f?Hv,r)j (r) 



(5.21) 

(5.22) 
(5.23) 



In the first of these equations we have used the Dirac equation to relate the 
bosonized wave functions /^o of the bound state to the fermionic components 
ho and j . 

As for the bound state contribution to the energy, this part of the moment 
of inertia is finite and will not be Pauli-Villars subtracted. Adding the c — c 
and b — c contributions the moment of inertia is given by 



M 2 



= 6 b - c + C " C (M) - — ^ C (A) . 



(5.24) 



Proceeding in an analogous way one can obtain expressions for the expec- 
tation values of other observables as well. Expressions for (E 3 ) and (L 3 ) in 
terms of mode sums have been derived in p^| . The spin expectation value is 
given by 

lN r ^ (lf>n\T3\lf>m)(ll>m\E3\li>n) 



(Es> 



ft 9 ' 
° z m,n 



E m E n 



(5.25) 
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It can be rewritten in terms of the Euclidean Green function as 



(E 3 ) = ~ f°° — I d 3 x I d 3 x'tr r 3 S^(x,x',u)a 3 Sl,(x',x,u) 



e 



in 



(5.26) 



where the contour again includes the bound state. The bound state and 
continuum contributions are obtained separately, as above. 



6 Observables 

The expressions for the energy and for the moment of inertia have already 
been presented. The mass of the low-lying baryons with K-spin is given by 
the sum of static and rotational energy as 

Mj = Eo + 4±jp±, (6.1) 

therefore 

M N = E + A (6.2) 

and 

15 3 3 

M A - M N = — - — = — . (6.3) 

89 89 29 y 1 

The nucleon sigma term is defined as 

E = m Q [ d 3 x(qq) . (6.4) 



As shown in [2{J it is given simply by the symmetry breaking part of the 
energy as 

E = E bl . (6.5) 



The experimental value is 45 ± 9 MeV |[30| . The pion- nucleon coupling con- 
stant can be obtained 0] from the long range behavior of the meson profile 

n t t 2 • ^ Q ^ exp(m^r) 

C — Jtt hm r sin » (6.6) 

r^oo 1 _)_ m7r? - 



as 



^nn = -7rM N C . (6.7) 
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The axial vector coupling constant is given [[H], by the expectation value 

9 a = (p T |7o 737375b T) • (6.8) 
It consists of a valence state contribution 



9a u = ~tr d 3 x( 7o7375 T 3 )Mx)4(x) 



3 

— - / drr 



h\r) - l -f{r) 



(6.9) 



and a continuum part which, using the Euclidean Green function, can be 
written as 

AT r roc rjii 

S/T = -^tr d 3 x ( l0l3l5T3 )S E (x,x,v) 

6 J J —00 Z71 

poo Ay poo 

■ 7T drr 2 Y,l(2K+l)s n (r,r,v)-(2K-l)s 22 (r,r,v) 
Jo Ztc Jo 

- (2K + 3) s 33 {r, r, v) + (2K + 1) s A4 {r, r, v) 



2 „ r r°° dv r°° 

— N r . 
9 



+A^K (K + 1) {s 23 (r, r, v) + s 32 (r, r, v)} . 

(6.10) 

Pauli-Villars subtraction is implied. With the Goldberger-Treiman relation 
the axial vector coupling constant can also be calculated via 

9*7* = ih^NN ■ (6.11) 

The quadratic radius of 

(R 2 ) hou = jT r 4 dr {h 2 (r) + j 2 (r)} , (6.12) 

1 rOO pOO 

(-R 2 )con = dv I r A dr V {s n + s 22 + s 33 + s 44 } (6.13) 

Z7T JO JO ~p 

has to be compared with an experimental value of 0.62 fm 2 . The continuum 
part is convergent, but so small that regularizing the integral does not change 
the result. 
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7 Numerics 



We have numerically implemented the expressions for the energy and its 
functional derivative presented in section 3 in the way described in || for the 
energy, and in |17| for its functional derivative. 

The iteration proceeds as follows: For a given meson profile $(r) one 
computes the mode functions and evaluates the functional derivative of the 
energy. One then requires the vanishing of the functional derivative, 



5E 



bou 



V 



+ 



6E hr M 2 5E ( 



A 2 




. 



(7.1) 



4>{Z)=Z${r) 



As can be seen from Eqs. (|3.23| ), (|3.27| ) and (|3.29|) , this equation takes the 
form 



5E 



<f)(Z)=Z0(r) 



z a [A(r) cos(i?(r)) + B(r) sin(tf(r))] 



(7.2) 



The coefficient functions A(r) and B(r) are the results of the numerical 
computation. Extremizing the energy by requiring the functional derivative 
to vanish fixes d(r) via tan(i?(r)) = —A(r)/B(r). This profile is used as 
the input for the next iteration. This method of iteration has been used 
previously by f3"2fl . 

The functions h a± (u,r) have been computed in order (1) and (2) by solv- 
ing (|4.4j) and its recursion (|4.7D using a Runge-Kutta scheme. The accuracy 
of these solutions was checked by using the Wronskian relation, which was 
constant to at least 6 significant units. The sum over the K-spin was ex- 
tended to -ft'max = 16 during the iteration and to K max = 20 for the final 
result. It is straightforward to derive, e.g., using asymptotic expansions for 
the lowest order perturbative contributions, that the power behaviour in an- 
gular momentum should be, after regularization, as K~ 3 . In Fig. 2a we show 
the terms in the sum over angular momenta for the integrand of the energy 
at v — 1 with a power fit AK~ 3 + BK~ A . An analogous example is given 
for the v integrand of g^, Eq. ( |6.10| ) at v — 0. Using this power fit it is 
then straightforward to include the sum above K max . So with a very good 
approximation the angular momentum sum runs, effectively, up to K = oo. 
Likewise the integral over v can be extended to v = oo by using a power fit. 
The v integrand for the energy is displayed in Fig. 2c. Here the expected 
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power behaviour is v~ 2 and, based on the power fit displayed in this figure 
we have appended the integral above v = 5M. 

The numerical results for the energy and other static parameters are 
presented in Table [I] and in Figs. 3-7. 



8 Results and conclusions 

We have presented here a self-consistent computation of the nucleon ground 
state in the Nambu-Jona-Lasinio model. In contrast to most previous cal- 
culations we have used a Pauli-Villars cutoff. It has been shown recently 
H that such a cutoff is favored by parton sum rules. As the main object 
of this work we have introduced an alternative method of numerical compu- 
tation, based on Euclidean Green functions instead of the use of the quark 
eigenfunctions for real energies. Our method has the advantage that it is 
not necessary to discretize the continuous spectrum of the sea quarks. While 
this discretization and the associated limiting procedure seem to be well un- 
der control, finite boundaries may introduce spurious effects, as discussed 
in ||27|| . Although our results essentially confirm those of other groups, this 
agreement is by no means guaranteed. 

Besides presenting the analytical framework for the computation of self- 
consistent profiles, based on explicit expressions for the energy and its func- 
tional derivative, we have also derived explicit expressions for other observ- 
ables. Again the quark sea contributions can be formulated in terms of the 
Euclidean Green function. 

Our numerical results are presented in Table [I] and plotted in Figs. 3 to 
7. In Table p] we also give some results obtained for g = 4 in Ref . || , when 
using the same regularization. In view of the difference of the numerical 
approaches the agreement is very satisfactory. This agreement holds as well 
for the mesonic profiles vr(r), plotted in Fig. 3. In Figs. 4 - 7 we plot the 
various parts of the energy, the axial vector coupling g^, the mesonic profiles 
i?(r) and the moment of inertia, as functions of the coupling g. One sees that 
a value of g ~ 4 is preferred by the comparison with experiment. 

The nucleon mass is still too high, but somewhat lower than the one 
obtained with Schwinger proper-time cutoff. A lower value for Mn can be 
obtained by minimizing the sum of fluctuation energy and rotational energy 
33| . We have not followed this issue here. It can be inferred from multiparti- 



cle dynamics (see, e.g. |34j) that one should subtract an energy corresponding 
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to the cm. motion of the quarks, which would result in a the mass near to the 
physical value. For the present selfconsistent condensate-quark state there 
is, however, no compelling proof for such a procedure. So such a subtraction 
is performed by some authors |35] , but has not become a general standard. 

A major difference to other publications on the chiral quark model is 
observed for the axial coupling constant g& which is computed from the quark 
current, and for which we obtain values between 1.15 and 1.37 as shown in 
Fig. 5. Most other authors have used the Schwinger proper-time cutoff, they 
obtain results well below 1. However, similar values for gA have been obtained 
recently by Golli et al. [2_SJ who use a Gaussian cutoff for the effective quark 
mass. The bound state contribution agrees with the one given, e.g., in ||. 
The major part of the increase of gA with g comes from the continuum 
part. The values for g\ obtained from the asymptotic behavior of the soliton 
profile via the Goldberger-Treiman relation are somewhat lower than those 
computed directly from the quark currents, they are in the range between 
1.15 and 1.17. While a small violation of the Goldberger-Treiman relation 
is expected, the appreciable increase of this violation with g is somewhat 
troublesome. 

As to the variation of the axial vector coupling with the coupling g, we 
find a monotonous increase. The bound state contribution decreases with g, 
this is, however, overcompensated by an increasing continum contribution. 
This trend is different from the one found by other authors (see, e.g., |5|, pSfl ), 
mostly using the Schwinger proper-time cutoff. It agrees, however, with the 
one found in || for their Pauli-Villars I cutoff, as we may infer from the two 
values gA = -94 at g = 3.85 and g^ = .96 at g = 4. The dependence of 
gA on g has again the same trend as ours in Ref. []22j with the Gaussian 
mass cutoff] In the same way the trend of the isoscalar quadratic radius (R 2 ) 
with g is found opposite to the one in other cutoff schemes; it again agrees 
with the Pauli-Villars I results of [|| . So it seems that it is the regularization 
scheme and not some numerical deficiency which causes the differences in the 
trends, an unsatisfactory situation which requires further investigation. 

While the absolute value of gA obtained here suggests a satisfactory agree- 
ment with experiment (see, e.g., [|36|) it has to be taken into consideration 
that gA is only calculated in O(Q ). As shown in the next orders in 



5 The correspondence between their parameters and ours is somewhat involved, in the 
range considered here their values of l/G are monotonuos with our coupling g, as are their 
and our values for the nucleon mass. 
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fl lead to additional contributions of g\ ~ 0.4 and g\ ~ 0.2, if computed 
with Schwinger proper-time regularisation. Thus in fact is overestimated 



here, as in nearly all regularization schemes, as discussed in |28[. The prob- 
lem is not resolved entirely, however, as apparently the expansion in Q does 
not converge well, so that higher corrections may still modify the results 
appreciably. 

In conclusion we have presented here a new approach to computing self- 
consistent meson profiles and static observables of the nucleon in the Nambu- 
Jona-Lasinio model, using a Pauli-Villars cutoff. The agreement with pre- 
vious analyses using the same cutoff but different numerical methods is sat- 
isfactory in general, some results given here are new. In view of the fact 
that our numerical procedure is rather economical we think that it is worth- 
while to pursue its application, e.g., to alternative versions of the model or 
to similar self-consistency problems. 
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A Partial waves in the K spin basis 

The expansion with respect to K-spin harmonics Ey |24j 




(A.l) 

written here for parity (— l)^ K+l \ reduces the Dirac equation to radial equa- 
tions for four coupled partial waves ttj. The Hamiltonian acting on the radial 
wave functions with parity (— via 

Hf Uj = Eu t (A.2) 
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is given by 



n 



K 1 



d 
dr 





K 



d 
dr 



K + 1 













r 



+ 



where 
















d K 










dr r 


C(r) 


—cS sS(r) 


° 1 


—cS(r) 


-C{r) 


-sS{r) 1 


sS(r) 


-C(r) - 


-cS{r) 





—sS{r) —cS(r) 


C{r) \ 


S(r) 


= Msin(tf(r)), 




Or) 


= Mcos($(r)), 





d 

dr 



K + 2 



(A.3) 



(A.4) 
(A.5) 

(A.6) 
(A.7) 



a = 2^K(K + 1)/(2K + 1), 
c = 1/(2^+1). 

We "square" the Dirac equation to obtain an effective Klein-Gordon equation 

f J = E 2 f l) (A.8) 



K 



-Ag + M 5ij + V l3 



where 



A 



K = S .±± r 2± _ K i(^ + 1) 

tJ lJ r 2 dr dr r 2 



(A.9) 



\K+1 



The orbital angular momenta K{ are given by K\ — K — 1, K2 = K3 = K 
and K4 = K + 1. The bosonized wave functions refer to the same K-spin 
basis (|A.1|) as the fermionic ones. The potential for the parity ( 
given by 

V K = 



is 



c S' + 



2KS\ 
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S'-- 
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-S' + 
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where 



S'(r) = M— sm(i?(r)) , 
dr 

C'ir) = m4- cos(0(r)). 
dr 



(A.ll) 
(A.12) 



For the parity (— 1) the sign of the mass has to be changed. 

While above we have formulated the Dirac and effective Klein-Gordon 
equations for wave functions with real minkowskian energies E we mainly 
work with euclidean mode functions whose argument we denote with v = 
—%E. So the mode equations are analogous to ( A.8 ) with E 2 — > — v 2 . In 



term of these mode fmctions the bosonic Green function is given by ( |3.21| ). 
The fermionic euclidean mode functions are related to the bosonized ones via 



U4 



{iv + C)h- ^A/ 2 -f 2 - cSf 2 + sSf 3 



{iv-C)f 2 



r 




K — 


1 


r 




K + 


2 


r 



fi + fi- cSfx - sSf 4 
fA-fi- cSfi + sSfi 



K 



(iv + C)U- -h + fs~ cSfc - sSf 2 



(A.13) 
(A. 14) 
(A.15) 
(A.16) 



for the parity (—lp K+1 \ thus the Green function becomes complex. With 
these radial functions the fermionic Green function Se(x-, x', v) in the K-Spin 
basis reads 

s nm (r, r', u)=k [9(r - /)<>, r)/*>, r') + 9(r' - r)<>, r)/£>, r') . 

(A.17) 



B Some relations for K-spin harmonics 

The action of r 3 , 03 and (X3T3 on the K-spin harmonics p| 



is 
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~K 
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~K 



z ^K,K Z 
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z ^K,K Z 
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A 2 - 
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K 2 
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(B.l) 
(B.2) 
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K x (2K-l)y/K " 2 (2K-1)K " 

jKj{K + iy-K^ K+1Kz _ K Z (2K-1) M 

(2K + 1)VK+1 ~ ^(2^ + 1) " 2 1 j 

, 9 \f R2 ~ K z W K-1,K Z 



(2K+1)JK(K + T) K(2K+1) 



^K+1,K Z . ^K,K Z 



+ 2 ; (B.7) 



(2K +1)(K+ lr 1 {2 K + 1)^K(K + 1) 
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2 (2K+1)(K + 1)- + (2K + l)J]^Tf 2 { } 



(2K + l)(K + ir 3 (2K + l)y/K{K + lf* 

((K + l)*-K*W + 2)*-K*) ^ 
(2K + 3)J(K + l)(K + 2) 



+2- 



K z ^(K + iy-K^ K+1Kz 
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Figure captions 



Fig. 1 : The complex z/-plane: Solid line: the contour C_ around the neg- 
ative continuum states; dashed line: the deformed contour along the real v 
axis, including, in addition, the bound state. 

Fig. 2a : Asymptotic behavior of partial wave sums and v integrands: 
we display the behaviour of the single contributions to the sum over partial 
waves for the v integrand of the energy at v — 1. The continous line is the 
power fit as described in the text. 

Fig. 2b : the same as Fig. 2a, partial wave contributions to the v inte- 
grand for gA at v — 0. 

Fig. 2c : the same as Fig. 2a, the v integrand for the energy. 

Fig. 3 : The mesonic profile of the nucleon for g — 4. We display the 
pion field n = f n sin(t?) as a function of r for g — 4.0 (solid line). For com- 
parison we also plot the result obtained previously by Doring et al. [5] using 
a different technique (dots). 

Fig. 4 : The energy of the nucleon. We display the energy as a func- 
tion of g (solid line), the bound state energy (dotted line), the zero point 
energy (dashed line) and the energy of the symmetry breaking term (long 
dashed line). For comparison we also plot the total energy and its parts 
obtained previously by Doring et al. [5] for g = 4 (dot, square, diamond and 
cross) . 

Fig. 5 : The axial vector coupling constant. We display the axial vector cou- 
pling constant as a function of g. The squares are obtained by calculating 
the expectation value, the dots by using the Goldberger-Treiman relation. 
The solid line is the experimental value. 

Fig. 6 : The mesonic profile for various couplings. We display the mesonic 
profile d as a function of r for different coupling constants, (g = 3.8 solid 
line, g = 4.0 dotted line, g = 4.2 dashed line, g = 4.4 long dashed line, 
g = 4.6 dash-dotted line). 
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Fig. 7 : The moment of inertia. We display the moment of inertia 9 as 
a function for different coupling constants (dots). The solid line is the ex- 
perimental value. 
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[MeV] 
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-Etotal 


[MeV] 


1157 


1145 


1131 


1116 


1099 


1139 








l/9 h - c 


[MeV] 


199 


218 


233 


247 


259 












i/e c - c 


[MeV] 


1404 


1264 


1141 


711 


575 












1/ #total 


[MeV] 


174 


186 


193 


183 


178 












M A - M N 


[MeV] 


261 


279 


290 


274 


267 






294 


M N 


[MeV] 


1222 


1214 


1203 


1184 


1166 






939 


(-R 2 )bou 
(-R )con 
(R 2 ) total 


[fm] 


0.56 


0.56 


0.58 


0.60 


0.62 












[fin] 


0.001 


0.002 


0.003 


0.004 


0.005 












[fin] 


0.56 


0.56 


0.58 


0.60 


0.62 


0.55 




0.62 


„bou 




0.72 


0.71 


0.71 


0.70 


0.70 


0.71 








„con 

9a 




0.43 


0.49 


0.55 


0.61 


0.67 












„total 

9a 




1.15 


1.20 


1.25 


1.31 


1.37 






1.23 


S'ttNN 




11.6 


11.7 


11.8 


11.8 


11.8 






13.1 | 
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G-T 

9a 




1.15 


1.16 


1.16 


1.16 


1.17 


0.96 


1.23 



Table 1: Static properties of the nucleon 
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